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Abstract 

The hyperfine induced 2s2p ^Pq, ^5*0 El transition probabilities of 

Be-like ions were calculated using grasp2K based on multi-configuration Dirac-Fock 
method and HFST packages. It was found that the hyperfine quenching rates are 
strongly affected by the interference for low-Z Be-like ions, especially for 2s2p ^Pq — > 
2s^ ^Sq transition. In particular, the trends of interference effects with atomic 
number Z in such two transitions are not monotone. The strongest interference 
effect occurs near Z = 7 for 2s2p ^Pq 2s^ ^Sq El transition, and near Z = 9 for 
2s2p 2s^ ^Sq El transition. 

PACS: 31.30.Gs; 31.15.V-; 31.15.ag. 

key words: Hyperfine induced transition; Hyperfine spectroscopy; Interference effect; Be-like 
ions. 

1 Introduction 

Hyperfine spectroscopy is very important tool in study of atomic and nuclear physics, 
especially which can be used to check fundamental interaction [H El El HI E] such as elec- 
tromagnetic and electroweak interaction with high accuracy, determine nuclear properties 
*Tel. No.: +86(0)931 7971148, Fax No.: +86(0)931 7971277, E-mail: Dongcz@nwnu.edu.cn 
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P El El El [in], develop atomic clock [III HI [131 [H [I5] , and so on. Recently one kind of 
hyperfine transitions, which is known as hyperfine induced transition or hyperfine quench- 
ing, attracts more attention owing to analyzing spectra [16] , determining isotopic ratios in 
stellar and diagnosing low-density plasma [HI HE] besides those reasons mentioned above. 

Many theoretical researches have been made to predict hyperfine induced transition 
probability based on diff'erent method. The most important problem in the calculations 
of hyperfine quenching rate is how to treat the hyperfine interaction and the interaction 
with electromagnetic field. In the past decades, three methods have been developed to 
solve such problem, i.e. perturbative method [20], complex matrix method [21], |22] 
and radiative damping method [23l [2l]. They have individual merit and shortage in 
actual physical problems as shown in Ref . [23l [21] . Here we used perturbative method 
to evaluate hyperfine quenching rate, that is, the computations of the probability follow 
the determinations of hyperfine structure. This approach is valid as long as radiative 
widths are smaller than the fine separation between concerned levels. The crucial point 
is to properly choice perturbative states involved in because the differences of results 
calculated by different approximate model are sometimes very large. 

For hyperfine induced 2s2p ^Pq 2s^ ^5*0 transition of Be-like ions there have been 
several studies. Marques et al. firstly computed the probabilities using their developing 
complex matrix method [23]. However, in their computational model they neglected an 
important contribution from the 2s2p ^Pi level, which lead to a relative large discrepancy 
compared to the later theoretical [T71 [21] and experimental value [iHl [26]. Once again, 
Brage et al. gave some of transition rates by means of perturbative method in order to de- 
termine the isotopic of composition and diagnose densities of low-density plasmas|17]. In 
their investigation the infiuence of 2s2p ^Pi level on hyperfine induced 2s2p ^Pq — ^ ^Sq 
transition of Be-like ions were indicated. But the calculations were restricted within those 
abundant elements in stellar. Later, Schippers et al. measured the probability for Be-like 
Ti ions using resonant electron-ion recombination method in the heavy-ion storage ring 
TSR of Max-Planck Institute for Nuclear Physics, Heidelberg, Germany [26j. The exper- 
imental result was almost 60% larger than Marques et al. theoretical value. This led to 
Cheng et al. renewedly calculating this decay rate along Be-like isoelectronic sequence by 
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perturbative and radiative damping method [21]. The latest theoretical results reduced 
the discrepancy to 20%. Their investigations shown again that the contribution from 
2s2p ^Pi transition amplitude and the interference effect from 2s2p 3Pi and 2s2p ^Pi 
perturbative states on hyperfine induced 2s2s 3Po ^5*0 transition of Be-like ions 

can not be neglected. Even though, it seems that characteristic about interference effects 
in hyperfine induced 2s2p ^Pq ^Sq transition of Be-like ions were still not very clear 

to be revealed. 

As we know, ^P2 level is another metastable state of sp configuration. It can decay 
to the excited state ^Pi by magnetic dipole (Ml) transition and to the ground state by 
magnetic quadurpole (M2) transition with large branch ratio. Many publications have 
been concerned with determination of the probabilities [271 EHl EHl [301 EI] to diagnose 
the low-density plasma and to probe relativistic and QED effects by accurate transition 
energy and probability. However, information about another transition process, that is, 
hyperfine induced El transition is very scarce. While, Gould et al. and Andersson et 
al. have pointed out significant influences of hyperfine induced El transition on the 
lifetime of ^P2 level in He-like |32j and Zn-like ions [33] , respectively. Dubau et al. have 
also shown that quantum interference between the El transition of hyperfine induced 
and M2 transition has obvious effects on increasing the degree of linear polarization of 
^P2 ^5*0 in He-like ions [H], which can affect the modelling and diagnostics of high- 
temperature astrophysical and laboratory plasmas with an anisotropic non-Maxwellian 
velocity distribution of energetic electrons. 

Based on these reasons mentioned above, we further investigated hyperfine induced 
2s2p ^P2, ^Pq 2.3^ ^Sq transition of Be-like ions in detail using grasp2K [3S] based 
on multi-configuration Dirac-Fock and HFST [36j packages. The regular of interference 
effects for hyperfine induced 2s2s ^Pqi 2s^ ^5*0 transition of Be-like ions was shown. 

In particular, it was found that the trends of interference effects with atomic number Z 
in such two transitions are not monotone. 
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2 Theory 



2.1 Wavefunction of hyperfine level 

In the present of hyperfine interactions, which couple the nuclear / and electronic 
J angular momenta to a total angular momenta F = I + J, only F and Mp are good 
quantum number other than the electronic angular momenta J and Mj. Then, the wave 
function for the system can be written by 

\FMF) = Y,h,jhJIFMF), (1) 

where h^^j is the mixing coefficients due to hyperfine interaction, and are obtained in first 
order perturbation theory as the ratio between the off-diagonal hj^perfine matrix elements 
and the unperturbed energy differences 



, _ {lJIFMF\Hnfs\lMFMF) 
EM) -EM ' 



the subscript labels the concerned level. The hyperfine interaction Hamiltonian Hhfs in 
this formula can be represented as a multipole expansion EH] , 

i^,^, = ^T^■M^ (3) 

fc>l 

where T^'') and M^^) are spherical tensor operators of rank k in the electronic and nuclear 
spaces, respectively [37]. In the following discussion we only include main the nuclear 
magnetic dipole (k=l) and electric quadrupole (k=2) interaction. Applying Racah's al- 
gebra, hyperfine interaction matrix elements can be further written by 



{jJIFMF\H^,f,\^oJoIFMF) ^ i-iy+''^'-^[i2J +l)i2I +1)]'^' { ^ ^ ^ (7J||T(i)||7Jo)(/||M(i)||/), (4) 

Jo / 1 

{^JIFMF\HhfsMIFMF) = {~ly+'''-^[{2J+l){2I+l)Y'^ ^ ^ W7J||t(2)||7Jo)(/||M(2)||/). (5) 

Jo / 2 

The reduce matrix elements of the tensor M*^*^) are related to the conventionally defined 
nuclear moment. 



(/||mW||/) = /.n/^, (6) 



(/||M(2)||/) = ^ 



(2/ + 3)(/ + l)_ 



2 ^ 7(2/ - 1 
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where is nuclear magnetic dipole moment in /Ua? of the nuclear magneton, and Q is 
electric quadrupole moment in barns. 



2.2 Hyperfine transition probability 

The electric dipole (El) transition probability between two different hyperfine levels 
\FMf) and \F'M'p) is given by 



A, ,3 

A=^i:\{FM,\q^^^\F'M'F)\\ (8) 
where Q'--'"-* is the electric dipole tensor operator. Substitute (1) into above formula, then 

A, ,3 1 

A = I E E KjK'J' ilJIFWQ^'^ \h'J'IF')\^ (9) 

7 J 7' J' 

Because operator Q^^^ only act on electronic parts, reduced matrix element {'yJIF\\Q^^^ \ \'j'J'IF') 
can be simplified based on Racah's algbra [39], 

(7J/F||Q«||yj7F') =(_i){J+/+F')^(2F + l)(2F'+l) ^ ^ M (7J||gW||yj'), 

^ [ F' J' 1 J 

(10) 

therefore, 

4, ,3 I J F / 1 

^ = ^(2F + l)|EEV/^7'^' (7^IIQ^'^ll7'/)P, 

(11) 

where u is the transition energy in Hartree. The reduced transition matrix elements of 
the electric dipole operator can be obtained as square roots of the corresponding line 
strengths. 

Using similar method one can obtain other type hyperfine induced transition proba- 
bility such as Ml, E2, etc. [HI [231 SOj. 

As can be seen from the derivation, hyperfine transition probability depends on nuclear 
parameters. It is not convenience for us to further discuss the trend of the rate along 
atomic number Z. Therefore, we generalized Brage et al. method [IT] so that hyperfine 



transition rate is independent of nuclear properties, which is called reduced hyperfine 
transition probability A^'. By defining reduced hyperfine mixing coefficient h'^\ 



h^i = L^i - (12) 

^4(l + /-l)(2/+l)]V2iy(/J0 J/;F1)' ^ ^ 



then, 



"^'^ 7^ I'J' \ F' J' 1 



(13) 



where W{IJo JF, Fl) are 6j-symbol in eq(4). To simplify we neglected electric quadrupole 
hyperfine interaction in above equation due to quite weak compared to the magnetic dipole 
interaction. 



2.3 Electronic wave function 

The electronic wave functions I7J) were computed using the grasp2K program package 
[35] . Here the wave function for a state labeled 7 J is approximated by an expansion over 
jj-coupled configuration state functions (CSFs) 

|7J) = EcA- (14) 
j 

The configuration state functions $j are anti-symmetrized linear combinations of products 
of Dirac orbitals. In the multi-configuration self-consistent field (SCF) procedure both the 
radial parts of the orbitals and the expansion coefficients are optimized to self-consistency. 
In the present work a Dirac-Coulomb Hamiltonian was used with the nucleus described 
by an extended Fermi charge distribution. 

Once the radial orbitals have been determined relativistic configuration interaction 
(CI) calculations can be performed. Here higher-order interactions may be included in 
the Hamiltonian. The most important of these is the Breit interactions 

_ _f;| °< ■ °.cosK.,) ^ ^^^^^^ ^^) CosK.,) - jj^j 

where photon frequency is obtained as the difference between the diagonal Lagrange 
multipliers and ej associated with the orbitals. However, this is invalid when shells are 
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multiply occupied, and the diagonal energy parameters of correlation orbitals with small 
occupation numbers may be large positive quantities totally unrelated to binding energies 
H2] . For this reason, the zero- frequency limit have been adopted in present calcu- 
lations. In the configuration interaction calculations the main quantum electrodynamics 
(QED) effects can also be included. 

Tensor algebra used for evaluating hyperfine and electric dipole matrix elements be- 
tween CI wave functions assumes that the wave functions are built from a common orbital 
set. This is a severe restriction since high-quality wave functions demands orbitals op- 
timized for the specific state. To relax this and to be able to compute matrix elements 
between wave functions built from independently optimized orbital sets, biorthogonal 
transformation techniques introduced by Malmqvist can be used jlSl HI] . 

3 Results and discussions 

3.1 Calculational model and method 

The accuracy of the calculated hyperfine induced transition rate depends on the 
number of perturbative states in Eq. (11) and on the accuracy of the electronic ma- 
trix elements. In practical calculation of hyperfine induced El transition probability for 
2s2p 3Po, ^Sq Eq. (11) is can be simplified to 

2 

(16) 

The differences for these two transitions concentrate on different hyperfine mixing coeffi- 
cient hs and transition energy u. As can be seen from Eq. (16), there exist interference 
effects caused by the ^Pi and ^Pi two transition amplitudes. 

The accuracy of the electronic matrix elements involved in above formula was deter- 
mined by wavefunction of electronic part. The main uncertainty in calculation of the 
wavefunctions comes from electron correlation effects, especially for low charged ions. For 
instance C^~*" ions, those matrix elements are fair sensitive to electronic correlation effects 
[lUj Wl\ 1^ . Therefore, it must be considered very carfully. As a starting point SCF 
calculations were done for the configurations in the even and odd state complex. In the 
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S=0,1 



calculations the wavefunctions of ls^2s^ ^5*0, ls^2s2p '^-Po,i,2 and ls'^2s2p ^Pi were deter- 
mined in extended optimal level (EOL) calculations [20], respectively. These calculation 
were followed by calculations with expansions including configuration state functions ob- 
tained by single(S) and double(D)-excitations from, respectively, the studied even and 
odd state reference configurations to active sets of orbitals n < 5. While for the 
active set was expanded to n < 7, and for N^"*" and to n < 6 in order to obtain 
satisfied results. The active sets were systematically increased allowing computed prop- 
erties to be monitored. Due to stability problems in the relativistic SCF procedure only 
the outermost layers of orbitals could be optimized each time. The SCF calculations 
were followed by CI calculations in which part of core-valence and core-core correlations, 
the frequency independent Breit interaction and QED effects was included. Finally, the 
hyperfine induced transition probability can be obtained by above mentioned theoretical 
method using HFST package [HS] . 

3.2 Hyperfine quenching rate 

In Table 1 we listed our calculated results for hyperfine induced transition probability 
as well as other theoretical and experimental value p!7l \2M l25l 126] . Nuclear parameters 
of concerned isotopes for Be-like ions were taken from ^T\. As can be seen from Table 
1, results of Marques et al. obviously deviate from others because two problems. One 
is that they neglected the contribution of ^Pi to hyperfine quenching rate of '^Po, and 
the other is that a ratio factor in transition rate was omitted in their method [23]. The 
present calculational results displayed in the forth column in Table 1 are in agreement 
with others for Z < 30, but not for ions with Z > 30. The reason for this is we neglected 
high order Breit interactions and QED effects in electronic wavefunction calculations, 
which lead to quite large errors in transition energy for high-Z ions. Due to limits of 
present method, we have not intrinsically solved this problem. Hence experimental [S2] 
and some of other accurate theoretical transition energies [21] were used to correct the 
hyperfine induced transition probability of ^Pq levle. Those corrected ones were presented 
in the fifth column of Table 1. It can be found that the consistence becomes obviously 
better with Brage et al. and Cheng et al. calculational value. While for Be-like ^°^Rh ion 



8 



the two order of magnitude difference is because different magnetic dipole moment was 
used. Based on this, we used the same method to correct other transition rates in next 
computations as welL 

In Table 2 we presented hyperfine quenching rate of in connection with corre- 
sponding transition energies used to correct those rates. Since angular momentum J of 
state unequal zero, it splits into several hyperfine levels, labeled by total angular 
momentum F. For individual hyperfine level satisfying select rule of electric dipole can 
occur hyperfine induced El transition, and this transition rate is dependent on angular 
number F and nuclear parameters. 

Due to lack of other results about hyperfine induced 2s2p 2s^ ^Sq transi- 

tion probability we can not make comparison. In order to confirm the rates is reliable, 
therefore, we further evaluated 2s2p 2s'^ ^Sq M2 and 2s2p 2s2j9 ^Pi Ml 

transition probabilities using our calculated line strength and experimental transition en- 
ergies. These results as well as other theoretical values were displayed in Table 3. As can 
be seen from this Table, the consistence among these results is quite good. It is indicated 
that our calculational hyperfine quenching rates of ^P2 level are credible. 

3.3 Interference effects in hyperfine quenching 

Brage et al. and Cheng et al. have dictated that interference effects strongly affect 
hyperfine quenching rate of ^Pq p!7l [2^ . From Eq.(16), we know that this effect occur in 
hyperfine quenching of ^P2 as well. To show clearly this phenomena, independent nuclear 
parameter of the transition amplitudes contributed from ^Pi and ^Pi were plotted in Fig. 
1. As can be seen from this picture, the interference effects work within a wide range of 
atomic number due to non-monotone change of transition amplitude contributed by ^Pi. 
It is interesting that the trend of transition amplitude with Z is similar between those 
two hyperfine quenching. An obvious difference is that the transition amplitude of ^Pi 
is dominant in 2s2p ^P2 — * 2s^ ^5*0 El transition while ^Pi for 2s2p ^Pq ^5*0 El 

transition. 

In order to reveal characteristic of the interference effects in such hyperfine induced 
2s2p ^Pq, ^P2 2s^ ^So transition of Be-like ions, it is convenience to define a function. 
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R , that proportion to the ratio between the two transition amphtudes. For example, for 
2s2p 2s^ ^Sq transition, 



According to the formula, the closer i?^' is to 0, the stronger interference effect is. There- 
fore, it clearly show the extent of interference effects. The trend of R^^ with Z for ^Pq 
and is plotted in Fig. 2, respectively. It worth noting from this picture that the 
interference effects in these two hyperfine quenching do not change monotonically and 
there exist minimum value for R^\^Pq) near Z — 1 and near Z = 9 for R^^i^P-i). Hence, 
the strongest interference effect occurs near Z = 7 for 2s2p ^Pq — > 2s^ ^5*0 El transition, 
and near Z = 9 for 2s2p — > 2s^ ^5*0 El transition. 

4 Conclusion 

In conclusion, we have calculated the hyperfine induced 2s2p ^Pq, ^P2 — > 2s^ ^5*0 El 
transition probability of Be-like ions using grasp2K based on multiconfiguration Dirac- 
Fock method and HFST packages. The interference effects resulted from ^Pi and ^Pi 
perturbative states in those two hyperfine quenching of 2s2p ^Pq, ^P2 2s^ ^Sq were 
studied in detail. It worth noting that the trends of interference effects with atomic 
number Z in such two transitions are not monotone. The strongest interference effect 
occurs near Z = 7 for 2s2p 3Po 2s'^ ^So El transition, and near Z — 9 for 2s2p P2 — > 
2s^ ^Sq El transition. 
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6 Figure caption 

Fig. 1. Transition amplitude of reduced hyperfine induced transition in a.u.. Left: 
2s2p transition; Right: 2s2p ^ Ss^ ^Sq transition. 

Fig. 2. The trend of i?^' for ^Pq and level with atomic number Z. Two blue arrows 
label the positions where inference effects are the strongest for the hyperfine induced 
2s2p ^Pq 2s^ ^Sq and 2s2p 2s^ ^Sq transition, respectively. 
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Table 1: Hyperfine induced 2s2p ^Pq 2s^ ^5*0 transition probability in and corresponding tran- 
sition energy in cm^^. The calculational results were compared with other theoretical and experimental 
value, uncorr. means the hyperfine induced transition probabilities was computed using present calcula- 
tional transition energy, while the corr. used the ones from NIST database [52] . 





Transition 


energy 


This work 










ions 


This work 


NIST [52] 


uncorr. 


corr. 


Ref.[T7] 


Ref.[ll] 


Ref.[55j 


Expt. 


13c 


52248 


52367 


8. 28 [-4] 


8.33[-4] 


9. 04 [-4] 


8. 223 [-4] 


2.00[-4] 




14N 


67251 


67209 


4.40 [-4] 


4.39 [-4] 


4.92[-4] 


4.40[-4] 


1.28 [-4] 


4[-4]±1.32'^ 


19p 


96666 


96590 


1.17[-1] 


1.17[-1] 




1.208[-1] 


3.60[-2] 






169054 


169802 


5.89[-2] 


5.97[-2] 


6.08[-2] 


6.011[-2] 


2.16[-2] 




39 Ar 


228716 


228674 


8.28[-l] 


8.27[-l] 










47Ti 


289562 


288190 


6.80[-l] 


6.71[-1] 




6.727[-l] 


3.56[-l] 


5.6[-l]^ 




352029 


348180 


4.98[-2] 


4.82[-2] 


5.45[-2] 


4.783[-2] 


3.27[-2] 




67Zn 


416600 


409827t 


5.00 


4.76 




4.732 


4.13 




85Rb 


537174 


523000 


43.3 


39.94 




39.35 


48.17 




103Rh 


693209 


661772t 


147.1 


128.0 




1.262 


1.91 




i3iXe 


903919 


843105t 


199.0 


161.5 




158.1 


262.67 





Cheng et al. [H] 
° Brage et al. [18] 
Schippers et al. [2S] 
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Table 2: Hyperfine induced 2s2p — > 2s^ ^Sq El transition probability A in s ^ associated with 



;ponding 


; reduced transition 


rate A 


in s 1 


and transition energy 


AE from 


NIST database 


ions 


AE 




F 


A 


ions 


AE 




F 


A 




52447 


9.87[-4] 


3/2 


7.30[-4] 




471780 


1.24 


3/2 


1.52[-1] 








5/2 











5/2 







67412 


2.52[-4] 


1 


2.08[-4] 


e^Zn 


640470 


4.40 


1/2 











2 


3.69 [-4] 








3/2 


8.81 








3 











5/2 


1.75[1] 


19-p 


97437 


1.21[-3] 


3/2 


1.26[-1] 








7/2 


1.79[1] 








5/2 











9/2 





2«Si 


177318 


1.70 [-2] 


3/2 


7.68[-2] 


85Rb 


1094800 


3.83[1] 


1/2 











5/2 











3/2 


1.85[2] 


39 Ar 


252683 


8. 28 [-2] 


3/2 











5/2 


3.65[2] 








5/2 


5.83[-l] 








7/2 


3.71[2] 








7/2 


1.03 








9/2 











9/2 


9.76[-l] 




2310547t 


4.08[2] 


3/2 


4.79[2] 








11/ 











5/2 







347420 


3.34[-l] 


1/2 





"iXe 


3785850 


4.96[3] 


1/2 


3.95[3] 








3/2 


4.82[-l] 








3/2 


1.30[4] 








5/2 


1.03 








3/2 


1.30[4] 








7/2 


1.17 


















9/2 
















^ Cheng et 
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Table 3: 2s 2p 2s^ ^ So M2 and 2s2p _> 2s2p ^Pi Ml transition probabilities of Be-like ions in s ^ connecting with corresponding transition energy from NIST 

database |52| in cm~^. The calculational results were compared with other theoretical and experimental value 



ions 




M2 { - 








Ml( 3p2 


- ''Pi) 






Transition enerj 


;y This work 


Ref. |30) 


Ref. [31] 


Transition energy 


This work 


Ref. [28] 


Ref. [29] 


Ref. PI 


13 c 


52447 


5.13[-3] 


5.190[-3] 


5.176[-3] 


56 


2.37[-6] 


2. 34 [-6] 




2.446[-6] 


14n 


67416 


1.14[-2] 


1.154[-2] 


1.147[-2] 


144 


4.03[-5] 


3.93[-5] 




4.070[-5] 


19p 


97437 


3.64[-2] 


3.678[-2] 


3.633[-2] 


587 


2.73[-3] 








28 Si 


177318 


2.42[-l] 


2.431[-1] 


2.410[-1] 


5174 


1.87 


1.83 






39 Ar 


252683 


7.89[-l] 


7.904[-l] 


7.858[-l] 


16820 


6.39[1] 


6.41[1] 


6.417[1] 




47Ti 


347240 


2.44 




2.4234 


42620 


1.03 [3] 




1.0369 [3] 




57Fe 


471780 


7.68 


7.652 


7.6459 


92655 


1.04[4] 


1.11[4] 






67Zn 


640470 


2.54[1] 




2.5289[1] 


180855 


7.52[4] 








85 Rb 


1094800 


2.28[2] 






480900 


1.31[6] 








103 Rh 


1996313t 


2. 90 [3] 




2. 8789 [3] 


1198404t 


1.86 [7] 








i3iXe 


3785850 


4.62[4] 






2758850 


2.10[8] 









t Cheng et al. [24] 
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